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ABSTRACT 

Aims. Redshift-space clustering anisotropies caused by cosmic peculiar velocities provide a powerful probe to test the gravity theory 
on large scales. However, to extract unbiased physical constraints, the clustering pattern has to be modelled accurately, taking into 
account the effects of non-linear dynamics at small scales, and properly describing the link between the selected cosmic tracers and 
the underlying dark matter field. 

Methods. We use a large hydrodynamic simulation to investigate how the systematic error on the linear growth rate, /, caused by 
model uncertainties, depends on sample selections and comoving scales. Specifically, we measure the redshift-space two-point corre¬ 
lation function of mock samples of galaxies, galaxy clusters and Active Galactic Nuclei, extracted from the Magneticum simulation, 
in the redshift range 0.2 < z < 2, and adopting different sample selections. We estimate /cg by modelling both the monopole and the 
full two-dimensional anisotropic clustering, using the dispersion model. 

Results. We find that the systematic error on /erg depends significantly on the range of scales considered for the fit. If the latter 
is kept fixed, the error depends on both redshift and sample selection, due to the scale-dependent impact of non-linearities, if not 
properly modelled. On the other hand, we show that it is possible to get almost unbiased constraints on /erg provided that the analysis 
is restricted to a proper range of scales, that depends non trivially on the properties of the sample. This can have a strong impact on 
multiple tracers analyses, and when combining catalogues selected at different redshifts. 

Key words, cosmology: theory, observations, dark matter, dark energy, large-scale structure of Universe 


1. Introduction 

The dynamics of the Universe is shaped by the gravity force. 
The understanding of gravity is thus the key to unveil the na¬ 
ture of the Universe and its cosmological evolution. After almost 
one century from its original formulation, Einstein’s General 
Relativity (GR) is still the dominant theory to describe gravity. 
However, despite its lasting experimental successes, some fun¬ 
damental open issues are motivating challenging efforts aimed at 
searching for possible e xpansions or radically alternative models 
dAmendola et al.ll2013h . Eirst of all, GR is not a quantum field 
theory and cannot be reconciled with the principles of quan¬ 
tum mechanics. Secondly, an increasingly large amount of as¬ 
tronomical observations cannot be explained by GR and bary- 
onic matter alone, requirin g the introducti on of a never observed 
form of dark matter (DM. IZwickvlll9.3^ . Finally, independent 
cosmological observations give indisputable evidences for an 
accelerated expansion of the Universe, requiring the introduc¬ 
tion of another dark component, generally du bbed dark energy 
(DE), when interprete d in the GR framework (iRiess et al.lll998t 
IPerlmutter et al.lll999ll . Whether GR provides a reliable descrip¬ 
tion of the large scale structure of the Universe, and thus the lat- 
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ter is dominated by unknown dark components, or, on the con¬ 
trary, such a gravity theory has to be corrected, or even aban¬ 
doned, is one of the key questions of modern physics and cos¬ 
mology. 


One of the most effective ways to test the gravity theory 
on large scales, that is where DM and DE arise, is to exploit 
the apparent anisotropies observ ed in galaxy maps, the so-calle d 
redshift-space distortions (RSD, lKaiseillT9^ lHamiltonlll998ll . 
Indeed, by modelling redshift-space clustering anisotropies it is 
possible to constrain the linear growth rate of cosmic structures, 
/ s d\og6/d\oga, where a is the dimensionless scale factor 
and 6 the linear fractional density perturbation, provided that 
the galaxy bias is known. The first RSD measurements were 
exploited primarily to estimate the mean matter density, Qm , 
that can b e derived directly fr o m f(z) when a gravi ty model is 
assumed jPeacoc k et al.l 2001; Hawkins et al.ll206^ . Later on, 
iGuzzo et al.l (12008 ) and Zhang et al.l (12008 ) showed that RSD 
can be effectively used to discriminate between DE and modified 
gravity scenarios. RSD thus started to be considered as a power¬ 
ful probe of gravity, and several applications to galaxy redshift 
surveys followed rapidly, both in the local Universe and a t larger 
redshifts, up to z ~ 1, su ch as 6dFGS (iBe utler et al.l2012h. SPSS 
( Samushia et al.l l2012t [C huang & Wang 2013 : [ch uang et ^ 
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VIPERS ( dgJa_ToiTe 6^1.1 120 1 3h . and BOSS dToieiro et alJ 
[20Tl[R^t alJl2012Ll20ll. 


All the above measurements have been obtained from large 
galaxy redshift surveys. First tentative studies are starting to con¬ 
sider also different tracers, but results are still strongly affected 
by statistical uncertainties, due to the paucity of the catalogues 
used. Thanks to the wealth of observational data now available, 
or expect ed in the next future, for instance from the ESA Euclid 
missiorQ dLaureiis et al.ll20flh . the NASA Wide Field Infrared 
Space Telescope (WFIRST) missior0 dSnergel et alJl201^ . and 
the extended Roentgen Surv^ with an Imaging Tel escope Array 
(eROSITA) satellite missiotO dMerloni et al.ll201^ . it will soon 
become possible to apply these methodologies on both cluster 
and Active Galactic Nuclei (AGN) catalogues. Galaxy clusters 
are the biggest collapsed structures in the Universe. They ap¬ 
pear more strongly biased than galaxies, that is their clustering 
signal is higher. Moreover, they are relatively unaffected by non¬ 
linear dynamics on small scales. This could help in modelling 
redshift-space distortions in their clustering pattern, as the effect 
of small scale incoherent motions, that generate the so-called 
fingers-of-God pattern, is substantially less severe relative to the 
galaxy clustering case, improving also the cosmological con- 
straints that can be extrac ted from baryon acoustic oscillations 
dVeropalumbo et al]l2014li . The main drawback is the low den¬ 
sity of galaxy cluster samples and the fact that these sources 
can be reliably detected only at relatively small redshifts. On 
the other hand, AGN can be detected up to very large distances. 
Powered by accreting supermassive black holes hosted at their 
centres, their extreme luminosities make them optimal tracers to 
investigate the largest sca les of the Universe, and thus its cosmo¬ 
logica l evolution (see e.g. lMarulli et al.l200^l200^lBonoli et alJ 
1 20091 and references therein). 

The aim of this paper is to investigate the accuracy of the 
so-called dispersion model (that will be described in |3]l as a 
function of comoving scales and sample selection. To achieve 
this goal, we exploit realistic mock samples of galaxies, galaxy 
clusters and AGN extracted from the Magneticum hydrodynamic 
simulation, at six snapshots i n the range 0.2 < z < 2. We adopt 
a similar methodology as in iBianchi et al.l (1201 2h . though sub¬ 
stantially extending that analysis. Specifically, i) in line with the 
majority of recent studies, we investigate the systematic errors 
of /cTg, instead of fi, ii) we consider different kinds of realistic 
mock tracers, instead of just DM haloes, hi) we analyse the de¬ 
pendency of the errors as a function of redshift, and iv) as a func¬ 
tion of different sample selections. In agreement with previous 
inves tigations (see e.g. IBianchi et al.l 120121 : iMohammad et alJ 
[2m^ , we find that the rough modelisation of non-linear dy¬ 
namics provided by the dispersion model introduces systematic 
errors on /erg measurements, that depend non trivially on the 
redshift and bias of the tracers. On the other hand, as we will 
demonstrate, it is possible to substantially reduce the systematic 
errors on /erg provided that the fit is restricted in a proper range 
of scales, that depends on sample selection. An alternative inves¬ 
tigation of the impact of the galaxy sample selection function on 
the ability to t est GR with RSD has been carried out recently by 
iHearinI (l2Q15h . analysing the small-scale effects of the assembly 
bias in velocity space. 

The software that we implemented to perform all the anal¬ 
yses presented in this paper can be freely downloaded at this 
link: http://apps.difa.unibo.it/files/people/federico.marulli3/ It 


' http://www.euclid-ec.org 
^ http://wfirst.gsfc.nasa.gov 
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consists of a set of C-I--I- libraries that can be used for sev¬ 
eral astronomical calculations, in particular to measure the two- 
point coiTelation fu nction and to model RSD (CosmoBolognaLih, 
iMarulli et al.ll201^ . We also provide the full documentation for 
these libraries, and some example codes that explain how to use 
them. 

The paper is organised as follows. In ^ we describe the 
Magneticum simulation used to construct mock samples of 
galaxies, clusters and AGN, whose main properties are reported 
in Appendix lAl In ^and ^we present the formalism used to 
model RSD and to measure the two-point correlation function in 
redshift-space mock catalogues, respectively. The results of our 
analyses are presented in ^ In ^we summarise our conclu¬ 
sions. Finally, in the appendices iBl and ICl we discuss the tests 
performed to investigate the robustness of the methods applied 
in our analysis. 


2. The Magneticum simulation 

The most direct way to investigate systematic errors on linear 
growth rate measurements from RSD is to exploit large numer¬ 
ical simulations, and to apply on them the same methodologies 
used for real data. As our goal is to study the effects of differ¬ 
ent sample selections, we make use of realistic mock catalogues 
of different cosmic tracers. Specifically, we analyse galaxy, clus¬ 
ter and AGN mock catalo gues extracted from the hyd rodynamic 
Magneticum simulatioifl (IDolag et ahllin preparation!) . 

This simulation is based on the paralle l cosmological 
TreePM-SPH code P-GADGET3 (ISpringelll2005h . The code uses 
an entropy-conser ving formulation of smoothe d-particle hydro¬ 
dynamics (SPH) (ISpringel & Hernauisll l2002h and follows the 
gas usin g a low-viscosity SPH scheme to properly track tur¬ 
bulence (IDolag et al.l l2005h . It also allows a treatment of ra¬ 
diative cooling, heating from a uniform time-dependent ultra¬ 
violet background and star formation with the associated feed¬ 
back processes. The latter is based on a sub-resolution model 
for the multiphase structure of the interstellar medium (ISM) 
(ISpringel & HernQuistjl200^ . 

Radiative cooling ra tes are computed by fo llowing the same 
procedure presented by IWiersma et al.l (l2009ll . We account for 
the presence of the cosmic microwaye background and of 
ultrayiolet/X-ray ba ckground radiation from quasars and galax¬ 
ies, as computed bv iHaardt & Mad^ (1200ill . The contributions 
to cooling from each one of the following 11 elements, H, He, C, 
N, O, Ne, Mg, Si, S, Ca, Fe, haye been pre-comp uted using the 
publi cly available CLOUDY photoionization code (IFerland et aO 
Il998h for an optically thin gas in photoionization equilibrium. 

In the multiphase model for star-formation 
(ISpringel & HernQuisti2003'9l . the ISM is treated as a two-phase 
medium where clouds of cold gas form from cooling of hot 
gas and are embedded in the hot gas phase assuming pressure 
equilibrium whenever gas particles are above a given threshold 
density. The hot gas within the multiphase model is heated by 
supernovae and can evaporate the cold clouds. A certain fraction 
of massive stars (10 per cent) is assumed to explode as super¬ 
novae type 11 (SNII). The released energy by SNII (10^' erg) 
is modelled to trigger galactic winds with a mass loading 
rate being proportional to the star formation rate to obtain a 
resulting wind velocity of Vwind = 350 km/s. Our simulation 
als o includes a det a iled m odel of chemical evolution according 
to iTornatore et al.l (l2007h . Metals are produced by SNII, by 
supernovae type la (SNIa) and by intermediate and low-mass 

^ www.magneticum.org 
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stars in the asymptotic giant branch (AGB). Metals and energy 
are released by stars of different mass by properly accounting 
for mass-dependent lifetim e s, with a lifetime function according 
to iPadovani & Matteuccil (Il993h. the m etallicity-dependent 
ste llar yields byiW ooslev & Weaveil ( 19951) for SNII, the yields 
by Ivan den Hoek & Groenew egenI ) 19971) for AGB stars and the 
yields bv iThielemann et al.ril2003l) for SNIa. Stars of different 
mass are initi ally distributed according to a Chabrier initial 
mass function (IChabrieiil2003l) . 

Most importantly, the Magneticum simulation also includes 
a prescription for black hole growth and for feedback from AGN. 
As for star formation, the accretion onto black holes and the as¬ 
sociated feedback adopts a sub-resolution model. Black holes 
are represented by collisionless “sink particles” that can grow 
in mass by accreting gas from their environments, or by merg¬ 
ing with other black holes. Our implemen tatio n is based on the 
model presented in ISnringel et al.l )l2005h and iDi Matteo et alJ 
(l2005h. including the same modifications as in the study of 
[Fabjane^L ( 20101) and som e new, minor changes, as described 


I Hirschmann et^ (I20l4) . 

The black hole (BH) gas accretion rate, Mbh, is es- 
timated by using the Bondi-Hoyle-Lyttleton approximat ion 

dHovle & Lvttletonll 19.391: iBondi & HovMI 1944lBondilll952h : 


iffsH - 


47rG iffgB./boostP 
(cj -I- 


( 1 ) 


where p and Cs are the density and the sound speed of the sur¬ 
rounding (ISM) gas, respectively, /boost is a boost factor for the 
density which is typically set to 100 and v is the velocity of the 
black hole relative to the surrounding gas. The black hole accre¬ 
tion is always limited to the Eddington rate. The radiated bolo- 
metric luminosity, Lboi, is related to the black hole accretion rate 
by 


Lboi - , 


( 2 ) 


where is the radiative efficiency, for which we adopt a 
fixed value of 0.1, standardly assumed for a radiatively effi¬ 
cient accr etion disk onto a non - rapidl y spinnin g black hole ac¬ 
cording to Shakura & Sunva^ (Il973h (see also Springelir2005t 
iDi Matteo et al.[2005). 

The simulation covers a cosmological volume with peri¬ 
odic boundary conditions initially occupied by an equal num¬ 
ber of 1526^ gas and DM particles, with relative masses that 
reflect the global baryon fraction, Qb/fiivi- The box side is 
896 /z“*Mpc. The cosmological model adopted is a spatially 
flat ACDM universe with matter density Dm = 0.272, baryon 
density Qb = 0.0456, power spectrum normalisation erg = 
0.809, and Hubble constant Hq - 70.4km s^'Mpe^/ cho¬ 
sen to match the sev en-year Wilkin s on M icrowave Anisotropy 
Probe (WMPA7, I Komatsu et al.l 1201 ll) . More details on 
the Magneticum simulation and th e deriv e d mock ca t alogue s 
can be found in Hirschmann et al.l ( |2014|) : ISaro et aP (l2014t) : 
iBocauet et al.l (1201 6l) . 

Firstly, we will analyse large samples of mock galaxies, clus¬ 
ters and AGN selected in stellar mass, cluster mass and BH mass, 
respectively. Then, we will consider several subsamples with dif¬ 
ferent selections, as reported in Tables lCniC.2l and lC.3l (see |5]). 


3. Redshift-space distortions 

To construct redshift-s pace mock cata l ogues , we a dopt the same 
methodology used in iBianchi et al.l (1201 2h and iMarulli et al.l 


(1201 iL l2012allbl) . We provide here a brief description. We con¬ 
sider a local virtual observer at z = 0, and place the centre of 
each snapshot of the Magneticum simulation at a comoving dis¬ 
tance, Dc, corresponding to its redshift, that is: 


Dc(z) = c 



(3) 


where c is the speed of light, and the Hubble expansion rate is: 


H(z) = Ho [Gm(1 + zf + Dk(l + zf 


+Qde exp 


I r 1 + Mz) 

I Jo 1+2 



(4) 


where Qk = 1 - Qm - w(z) is the DE equation of state, and 
the contribution of radiation is assumed negligible. In our case: 
Dk = 1, w{z) = 1, so Gde = = 1 - and the Hubble 

expansion rate reduces to: 


H(z) = Ho [Qm( 1 + z)^ + Da] . 


(5) 


To estimate the distribution of mock sources in redshift- 
space, we then transform the comoving coordinates of each ob¬ 
ject into angular positions and observed redshifts. The latter are 
computed as follows: 


V|| 

Zobs = Zc + —(1 + Zc) , (6) 

c 

where Zc is the cosmological redshift due to the Hubble reces¬ 
sion velocity at the comoving distance of the object and vy is 
the line-of-sight component of its centre of mass velocity. In this 
analysis we do not include either errors in the observed redshift, 
or geometric distortions caused by an incorrect assu mption of 
the background cosmology (see IMarulli et al.ll2012bl for more 
details). 

In the linear regime, the velocity field can be determined 
from the density field, and the amplitude of RSD is proportional 
to the parameter/, defined as follows: 

fiz) 

/(z) = ^ , (7) 


where the linear growth rate, /, can be approximated in most 
cosmological frameworks as: 

fiz) - Qj;,(z), (8) 

with y ^ 0.545 in ACDM. The linear bias factor can be esti¬ 
mated as: 


biz)^{- 


I /(gz) 

I fominz) 




(9) 


where the mean is computed at sufficiently large scales at which 
non-linear effects can be neglected. 

In the linear regime, the redshift-space two-point correlation 
function can be written as follows: 

f’'(s,p) = ^o(s)Po(m) + f 2 (s)P 2 (M) + ^4(s)P4(m) , (10) 


where p = cos 6 = jy/j is the cosine of the angle between the 


separation vector and the line of sight, s - ^sj_ + and P; 
are the Legendre polynomials (lKais eiill987t Lilie & Efstathioiil 


ll989l:lMcffllllll99()l:lHamiltonll 992HFisher et al.ll994 . Eq. (fTOl) 
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is derived in the distant-observer approximation, that is reason¬ 
able at the scales considered in this analysis. The multipoles of 
^(ji,i||)are: 




(bo-sf + -fa-% ■ bos + -^(fos) 




^dmW 


f/T. 


bo-g + -{fosY 


^DM(r) ^dm(^) 


f4(.) = 

8 




35 


(fcrgY 


^DM(r) 


crj 


££pmW _ 7 ^dm(0 
2 


crj 


o-t 


(lla) 

(llb) 

(12a) 

(12b) 

(13a) 

(13b) 


^ dm ('‘) = 4 r dr'^ouir'y 
Jo 

^dm('‘) = 4 r '^'■'^DM(r')'‘' 
Jo 


,'2 


,'4 


^ll) = J" dvfivy S|| - 


v(l +z) ) 

/tin 


/(v)- 


1 


cr 12 


V2 


exp 


V 2 |v| 
Cr 12 


dPavis & PeeblesI 1 19831 : iFisher et ^ 11994 IZurek et ^ 1 1994 . 
The quantity an can be interpreted as the dispersion in the pair¬ 
wise random peculiar velocities, and is assumed to be indepen¬ 
dent of pair separations. 

The dispersion model given by Eqs. llT0l) - (fT7l) depends on 
three free quantities, fag, bag and cri 2 (since ^dm af), and on 
the reference background cosmology used both to convert angles 
and redshifts into distances and to estimate the real-space DM 
two-point correlation function. This model has been widely used 
in the past years to estimate the linear growth rate both in config- 


Peacock et al.l 20011: Hawkins et alJ 1200.31: iGuzzo et al.l 

2QQ81 

Ross et al.l 

20071 

Cabre & Gaztanagal l2009a Ibl: IContreras et al. 

2013) and in Fourier space, by modelling the power 

spec- 


where ^(r) and fDM(r) are the real-space undistorted correlation 
functions of tracers and DM, respectively, whereas the barred 
functions are: 


(14) 


(15) 


(16) 


where the pairwise velocity v is expressed in physical coordi¬ 
nates. 

Since we are not co nsidering redshift err ors, we use the ex¬ 
ponential form for /(v) (IMamlli et al]l2012hh . namely: 


(17) 


trum (e.g. lPercival et al.]l2004lTegmark et alJl2004 iBlake et al.l 
1201 ih . In this analysis, we investigate the accuracy of the disper¬ 
sion model in configuration space, while alternative RSD models 
will be analysed in an upcoming paper. 

4. Methodology 

The aim of this work is to test widely-used statistical methodos 
to model RSD and to extract constraints on the linear growth 
rate. Thus, we will use methodologies that can be applied di¬ 
rectly to real data. On the other hand, we will not consider any 
specific mock catalogue. In order to keep our analysis general, 
that is not restricted to any specific real survey. To measure the 
two-point correlation func t ions o f our mock samples we make 
use of the lLandv & Szala"^ (Il993h estimator: 




00(r) - 20R(r) + RR(r) 
RR{r) 


(18) 


Eqs. (Illbl) . (I12bl) and (I13bl) are derived from Eqs. (II lab . (I12all 
and (I13al) respectively, using Eq. (O. In this analysis, ^DM(r) 
is estimated by Eourier transforming the linear power spect rum 
computed with the software CAMS (iLewis & Bridlell2002h for 
the cosmological model here considered. We could alternatively 
measure the DM correlation function directly from the snapshots 
of the simulation. However, by using CAMB we can substan¬ 
tially reduce the computational time, we get a smooth model 
with no errors due to measurements and, more importantly, we 
can closely mimic the analysis we would have done with real 
data. 

Eq. ([lai is an effective description of RSD only at 
larg e scales, where non-linear effects are negligible (but see 
e.g.lSco ccimarrol |2Q04 lTaruva et al]l201(]l:[Seliak & McDonald! 
1201 ll: I Wang et akboidl for a more accurate modelling). An em¬ 
pirical model that can account for both lin ear and non-linear 
dynamics is the so-ca l led dispersion model ([Peacock & DoddsI 
1 1994 IPeeblesI 119801: iDavis & PeeblesI 1 198^ . that describes 
the redshift-space correlation function as a convolution of the 
linearly-distorted correlation with the distribution function of 
pairwise velocities, /(v): 


where 00(r), OR(r) and RR(r) are the fractions of object-object, 
object-random and random-random pairs, with spatial separa¬ 
tion r, in the range [r - drjl, r + i5r/2], where 6r is the bin size. 
The random catalogues are constructed to be three times larger 
than the associated mock samples. As we are considering mock 
catalogues in cubic boxes, with no geometrical selection effects 
and periodic conditions, we could easily estimate the two-point 
correlation function directly from the density field, or comput¬ 
ing the random counts analyti cally. Nevertheless, we prefer to 
use the lLandv & Szala"^ (Il993ll estimator to mimic more closely 
the analysis on real data, as stressed before. In any case, this 
choice does not impact the results of our analysis. Moreover, as 
we verified, the number of random objects used in this work is 
large enough to have no significant effects on our measurements 
(see Appendix |B]l. We compute the correlation functions up to 
a maximum scale of r = 50 /r ' Mpc, both in the parallel and 
perpendicular directions. As we have verified with a sub-set of 
mocks, going to higher scales does not change our results, as the 
constraints on fag are mainly determined by the RSD signal at 
smaller scales. 

To model RSD and extract constraints on fag, we exploit 
the model given by Eqs. (IT0l) - (fT7]) . We consider two case stud¬ 
ies. The first one consists in modelling only the monopole of 
the redshift-space two-point correlation function at large scales, 
via Eq. (II Ibb . that is in the Kaiser limit, assuming that the non¬ 
linear effects can be neglected at these scales. In this case, the 
bag factor has to be fixed a priori. When analysing real cata¬ 
logues, this factor can be dete rmined either directly with the de¬ 
projection technique (see e.g. iMarulli et ak 2^^^j_orfrom_the 
three-point correlation function (see e.g. iMoresco et al.ll20l4 . 
or by combining clustering and lensing measurements (see e.g. 
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ISereno et alJl2015h . To minimise the uncertainties and systemat- 
ics possibly present in the above methods, we derive bcr% directly 
from the real-space two-point correlation functions of our mock 
tracers, through Eq. (O, while the DM correlation is obtained by 
Fourier transforming the linear CAMS power spectrum, fixing 
CTg at the value of the Magneticum simulation. 

The second method considered in this work consists in fully 
exploiting the redshift-space two-dimensional anisotropic corre¬ 
lation function ^n). An alternative approach would be to 
use the multipole moments of the correlation function. The ad¬ 
vantage is to reduce the number of observable s, thus helping in 
the computation of the co variance matrix fe.g. Ide la Torre et aH 
l2013tl^tracca et al.l2016ll . However, we are not mimicking here 
any real redshift survey. Thus, a detailed treatment of the statis¬ 
tical errors is not necessary for the purposes of the present pa¬ 
per, that focuses instead on the systematic uncertainties caused 
by a poor modelisation of non-linear effects. Therefore, we will 
use just diagonal covariance matrices, whose elements are es- 
timated analytically from Poisson statistics. As discussed in 
iBianchi et al.l (l2012b . the effect on growth rate measurements 
caused by assuming negligible off-diagonal elements in the co- 
variance matrix is small (of the order of only a few percent in 
that analysis), thanks to the large volumes of the mock samples 
with respect to the scales used for the parameter estimations. 
Appendix|C]provides some further tests, using the mock samples 
analysed in this wor k, that appear in over all agreement with th e 
iBianchi et al.l (12012h conclusions (see also lPetracca et al.ll2016h . 
Moreover, to investigate the accuracy of the dispersion model as 
a function of scales, it is more convenient to exploit f(5x, ■sy) ■ 
Indeed, following this approach, we can explore the full ry - rj^ 
plane, searching for the optimal region to be used to minimise 
model uncertainties. 

5. Results 

In this section, we present our main results obtained from mock 
samples of three different tracers - galaxies, clusters and AGN 
- using the two approaches described in ^ For each cata¬ 
logue, we analyse six snapshots corresponding to the redshifts 
z - {0.2,0.52,0.72,1,1.5,2). Moreover, at each redshift we con¬ 
sider different sample selections. In total, we analyse 270 mock 
catalogues, whose main properties, including the number of ob¬ 
jects in each sample, are reported in Tables lCniC.2l and lC.31 and 
in Fig.[T2]in Appendix lAl 

5. 1. The linear growth rate from the clustering monopole 

We start analysing the spherically averaged two-point correla¬ 
tion function - the clustering monopole - at large linear scales. 
The aim of this exercise is to investigate the accuracy of the RSD 
model in the simplest case possible, that is in the so-called Kaiser 
limit. 

Fig. [1] shows the ratio between the redshift-space and real- 
space two-point correlation functions, ^{s)l^(r). Specifically, we 
show here the results obtained at z = 0.2, for galaxies with 
log(MsTAR[/* * Mq]) > 10, clusters with log(M5oo[/!^* Mq]) > 
13, and AGN with Mq]) > 6.3, as reported by the 

labels. These correspond to the upper left samples reported in 
Tables ICniC.2l and 1C. 31 Results obtained from the other mock 
samples are similar, but more scattered due to the lower den¬ 
sities. The redshift-space correlation functions have been mea¬ 
sured directly from the simulation, through the procedure de¬ 
scribed in ^ The real-space correlation functions have been es¬ 
timated in two different ways. Either they are measured directly 
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Fig. 1. The ratio between the redshift-space and real-space two- 
point spherically averaged correlation functions, ^(s)l^(r), at 
z = 0.2, for galaxies with log(MsTAR[/!^^ Mq]) > 10 (top panel), 
clusters with log(M 5 oo[/!^* Mq]) > 13 (centralpanel), and AGN 
with log(MBH[/j^* Mq]) > 6.3 (bottom panel). The redshift- 
space correlation functions have been measured directly from 
the simulation, while the real-space ones have been either mea¬ 
sured from the simulation (red dots), or derived from the CAMB 
power spectrum (blue squares ). The error bars represent the sta¬ 
tistical noise as prescribed by Mo et al.l (Il992h . The black solid 
lines show the expected values of ^(s)l^(r) at large scales pre¬ 
dicted by Eq. dl Ibk The dotted blue and dashed red lines are 
the best-fit ratios estimated from the blue and red points, respec¬ 
tively. 


from the simulation, or they are derived from the linear CAMB 
power spectrum, and assuming a linear bias factor. The differ¬ 
ences at small scales, r < 5 Mpc, between the ratios com¬ 
puted with the measured (blue squares) and CAMB (red dots) 
real-space correlation functions are due to non-linear effects. 
We verified fhaf using th e non-linear CAM B power spectrum, 
via the HALOFIT routine dSmith et alJl200.3h . does not fully re¬ 
move the discrepancy. Nevertheless, as we model here the large 
scale clustering, this has no effects on our results. On the other 
hand, the small discrepancies at scales r > 40 /z ' Mpc can in¬ 
troduce systematics, that however result smaller than the esti¬ 
mated uncertai nties. The error bars in Fig. [1] show the statistical 
Poisson noise (iMo et alJll992h . Scales larger than 60 /z * Mpc, 
not shown in the plot, are too noisy to affect the fit. More specif¬ 
ically, we verified that a convenient scale range to get robust 
results is 10 < r[h^^ Mpc] < 50. The black solid lines show 
the expected values of ^(s)l^(r) at large scales, as predicted by 
Eq. (fTTbll . 

By modelling the clustering ratio, ^(s)l^(r), via Eq. dllall . 
it is possible to estimate the factor [3. We do this by measur¬ 
ing both ^(s) and ^(r) from the simulation. The result is shown 
by the dotted blue lines in Fig. [1] Instead, to estimate /erg 
from the monopole of the two-point correlation function, we use 
Eq. (II Ibl) . In this case the factor ber^ has to be fixed. This re- 
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Fig. 2. Top panel: the best-fit values of f(z)(T^iz) for galax¬ 
ies with log(MsTAR[/t * Mq]) > 10 (red dots), clusters with 
log(M 5 oo[/!^* Mq]) > 13 (blue squares), and AGN with 

log(MBH[/r^' Mq]) > 6.3 (green diamonds), obtained with the 
method described in ^5.11 The black line shows the function 
'CrgCz), where Qm( 2) and cr8(z) are the known values of 
the simulation. Bottom panel: the percentage systematic errors 

on f(z)0-s(z), [(/cr8)™"™''‘=A(/o-g)simulation]/(y^gpimulation . 

that is the percentage differences between the points and the 
black line shown in the top panel. The error bars have been es¬ 
timated with Eq. ( fT9b . The values reported have been slightly 
shifted for visual clarity. To guide the eyes, the light and grey 
shaded areas highlight the 5% and 10% error regions, respec¬ 
tively. 


quires to estimate the real-space two-point correlation function 
of the DM. As described above, we obtain the latter by Fourier 
transforming the CAMS power spectrum. Moreover, to estimate 
the mean linear bias, a range of scales has to be chosen in Eq. ([9|. 
We compute the linear bias in the range r < r[h-^ Mpc ] < 50, 
where r is a free parameter that sets the minimum scale be¬ 
yond which the bias is fairly scale-independent. Such a mini¬ 
mum scale depends on redshift and sample selection, and we 
estimate it for each sample considered. The dashed red lines are 
the best-fit ratios estimated with this method. 

As it can be seen by comparing red and blue lines, the two 
methods provide concordant results, both of them in agreement 
with the expectations. Indeed, the real-space correlation function 
provided by CAMB and assuming a linear bias is in reasonable 
agreement with the one measured directly from the simulation. 

The constraints on fcr^ as a function of redshift are shown in 
the upper panel of Fig. |2 As in Fig. [1] we show the results ob¬ 
tained for galaxies with log(MsTAR[/i ^ M©]) > 10, clusters with 
log(M5oo[/!^' Mq]) > 13, and AGN with log(MBH[/!^* M©]) > 
6.3. The black line shows the function ■ cr8(z), where 

Qm(z) and cr8(z) are the values of the Magneticum simulation. 
In the lower panel we show the percentage systematic errors on 
/erg, defined as [(/o-8)““™'"‘‘-(/cr8)“™‘"‘‘™]/(/cr8)“™‘"‘‘°'‘-100. 
The error bars have been esfimafed by propagating on /erg the / 


Fig. 3. The iso-correlation contours of the redshift-space 
two-point correlation function, corresponding to the val¬ 
ues [0.05,0.1,0.2,0.4,1,3], for galaxies with 

log(MsTAR[/i~* M©]) > 10, at z = 0.2 (black contours). The 
dot-dashed green and solid red contours show the best-fit model 
given by Eq. (UHl, with the real-space correlation function ^(r) 
measured from the simulation, and estimated from the CAMB 
power spectrum, respectively. The blue dashed contours show 
the linear best-fit model given by Eq. (ITOll with the CAMB real- 
space correlation function. 


errors provided by the scaling formula presented in iBianchi et al.l 
(l2012h . that gives the statistical errors as a function of bias, b, 
volume, V, and density, n: 

„ 9 ) 

where no ^ 1.7- IQ-^^ Mpe^^ and C = 4.9 ■ 10^ Mpc* ^ In 
this case, these error bars have to be considered just as lower lim¬ 
its, as the scaling formula has been calibrated based on fits of the 
full two-dimensional anisotropic correlation function iy). 
Moreover, they have been computed in different scale ranges 
and with Friends-of-Friends DM haloes, differently from this 
analysis, where we consider tracers hosted in DM sub-haloes. 
Using the full covariance matrix in this statistical analysis does 
not change our conclusions, as shown in Appendix IQ 

As demonstrated by Figs.[T]and|2l it is indeed possible to get 
almost unbiased constraints on ferg from the monopole of the 
two-point correlation function, at all redshifts considered, and 
for both galaxies, clusters and AGN, provided that the linear bias 
is estimated at sufficiently large scales. The advantage of this 
method is the minimal number of free parameters necessary for 
the modelisation. However, the linear bias has to be assumed, or 
measured from other probes. 

5.2. The linear growth rate from the anisotropic 2D ciustering 

To extract the full Information from the redshift-space two-point 
correlation function, the anisotropic correlation iy) or, al¬ 
ternatively, all the relevant multipoles have to be modelled. As 
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Fig. 4. Top panel of the upper window: the best-fit values of 
fiz)crg{z) for galaxies with log(MsTAR[/!”* Mq]) > 10. The open 
and solid blue squares show the values obtained by fitting the 
data with the models given by Eqs. Col) and d, respectively, 
in the scale ranges 3 < Mpc] < 35 and 3 < Mpc] < 

35. The open and solid red dots have been obtained with the 
method described in ^5.2.11 fitting the data with the models 
given by Eqs. (fTOl) and (IThl) . respectively. The black line shows 
the function ■ crg(z), where Qm(z) and o-s(z) are the 

known values of the simulation. Eor comparison, the grey dots 
show a set of recent observational measurements from galaxy 
surveys (see ^5.2.11) . Bottom panel of the upper window: the per¬ 
centage systematic errors on fiz)o-g(z), defined as 

■ 100, that is the percentage differ¬ 
ences between the points and the black line shown in the top 
panel. To guide the eyes, the light and grey shaded areas high¬ 
light the 5% and 10% error regions. Lower window: best-fit val¬ 
ues of r™" and of the method described in ^5.2.11 The grey 
areas show the regions explored by the method. 


discussed in ^ we consider the first approach. Differently from 
the analysis of ^5.11 here we can jointly constrain the two terms 
/cTg and bcTg. In this section we investigate the accuracy of the 
constraints on /erg provided by the dispersion model, as a func¬ 
tion of sample selection, with berg and cri 2 as free parameters. 


Fig. 5. The percentage systematic errors on f {z)crg{z) for galax¬ 
ies, [(/crg)™“™‘'=d-(/o-8)“™i^‘i°"]/(/crg)“™‘=‘“™ ■ 100, as a func¬ 
tion of redshift and sample selections, as indicated by the labels. 
S1 - S 5 refer to the five selection thresholds reported in Tables 
IC.lllC.2l and lC.31 for the different tracers and properties used for 
the selection. 


5.2.1. Galaxies 

We start presenting our results for the galaxy mock samples. 
The black lines of Eig. [3] show the iso-correlation contours of 
the redshift-space two-point correlation function, /(ij., ^n), for 
galaxies with log(MsTAR[/* * Mq]) > 10, at z = 0.2. The other 
lines are the best-fit models obtained with three different meth¬ 
ods. The dot-dashed green lines are obtained by using the full 
non-linear dispersion model given by Eq. d), with the real- 
space correlation function f{r) measured directly from the simu¬ 
lation. The blue dashed and red solid contours are obtained with 
/(r) estimated from the CAMS power spectrum, and by fitting 
the linear model given by Eq. (ITOl) and the non-linear one by 
Eq. d, respectively. The fitting is done in the scale ranges 
3 < Mpc] < 35 and 3 < /■||[/z“' Mpc] < 35. 

As one can see, while all the three models provide a good 
description of the anisotropic clustering at scales larger than ~ 5 
' Mpc, none of them can reproduce accurately the fingers-of- 
God pattern at small scales. This is expected for the blue con¬ 
tours, obtained without modelling the non-linear dynamics. For 
the other two cases, the discrepancy is due to the not sufficiently 
accurate description of non-linearities provided by the disper¬ 
sion model. We verified that the minimum scales considered for 
the fit have a negligible impact on the fingers-of-God shape, that 
is we find the same results even considering scales smaller that 
3 Mpc in the fit. The green contours appear in better agree¬ 
ment with the measurements at small scales, due to the fact that 
the real-space correlation function is measured from the mock 
cataloguqj. 


^ The model corresponding to the green contours provides constraints 
on not on /cg. We show it here only to highlight the differences at 
small scales between the dispersion model estimated with ^(r) measured 
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Fig. 6. The iso-correlation contours of the redshift-space 
two-point correlation function, corresponding to the values 
■Sii) = [0-2,0.4,1,3], for clusters with log(M 5 oo[/!^' Mq]) > 
13, at z = 0.2 (black contours). The other contours are as in 
Fig.E 


The best-fit values of /erg as a function of redshift are shown 
in the top panel of the upper window of Fig. |4l The open 
and solid blue squares show the values obtained by fitting the 
data with the models given by Eqs. (flOl) and (fT^ . respectively, 
that is they correspond to the blue and red contours of Fig. [3] 
(at z - 0.2). The bottom panel of the upper window shows 
the percentage systematic errors on /erg, that is [(/crg)™‘=‘‘®“‘^‘^‘^- 
(/cr8)^‘™''“‘“]/(/crg)'*‘™“™ ■ 100. 

As it can be seen, the best-fit values of /erg result strongly 
biased with respect to the true ones, with a systematic error 
that depends on the redshift. At z = 0.2, we get an underesti¬ 
mation of ~ 10 (20) % with the non-linear (linear) dispersion 
model, at z = 1 the error reduces to ~ 5 (8)%, while at z = 2 
we get an overestimation of ~ 20%, with both linear and non¬ 
linear m odelling. This result is fairly in agreement with what 
found by iBianchi et al.l (l2012h at z = 1, although the two anal¬ 
yses are not directly comparable, as IBianchi et al.l (1201 2l) con¬ 
sidered a sample of Friends-of-Friends DM haloes, slightly af¬ 
fected by fingers-of-God. As expected, the convolution given 
by Eq. (fT^ reduces the systematic error on /erg, especially at 
low redshifts. However, it does not totally rem ove the discrepan¬ 
cies, i n agreement with previous findings (e.g. [Mohammad et al.l 
1201 6l and references therein). Eor comparison, the grey dots 
of Fig. 0] show a set of recent observational measurements 
of fcTg from large galaxy surveys: 6dFGS at z = 0.067 
dBeutler et alj201^: SDSS(DR7 ) Luminous Red Galaxies at z = 
0.25,0.37 ( Samushia et alJl2012h from scales lower than 60 and 
200 /r 'Mpc; B OSS at z = 0.3,0.57,0.6 (iToieiro et al.ll2012l: 
Reid et aP 2012h : Wiggle Z at z = 0.44,0.6,0.7 3 ( Blake et^ 


2012 ): VIPERS at z = 0.8 (Ide la Torre et al.ll2013h . These results 


have been obtained using different RSD models, most of them 


from the simulation, and the one with ^(r) from the linear CAMB power 
spectrum, assuming a linear scale-independent bias. 
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Fig. 7. The best-fit values of /(z)cr8(z) for clusters with 
log(M 5 oo[/!^* Mq]) > 13. All the symbols are as in Fig. 01 


more accurate than the dispersion model considered in this work. 
Nevertheless, many of these measurements un derestimate ferg 
with r espect to GR-i-ACDM predictions (see e.g. lMacaulav et aO 
l2013h . In line with our findings, this could be explained, at least 
partially, by model uncertainties still present in the more sophis¬ 
ticated RSD models considered. 


The direct way to reduce these systematics is t o improve the 


modelisation of RSD at non-linear scales ("see e.g.|Scoccimarro 

2004 

:lTaruva et al.ll2010l: Seliak & McDonaldll201 It Wan? et al. 

2014 

; Ide la Torre et al.||2013|). We explore here a different ap- 


proach, investigating the dependency of the systematic error on 
the comoving scales considered in the analysis. In the forth¬ 
coming plots, we show the results obtained by repeating our 
fitting procedure for different values of the minimum perpen¬ 
dicular separation and the maximum parallel separation used in 
the fit, that is r™" and r™“. By increasing the value of we 
can cut the region more affected by fingers-of-God anisotropies, 
while by reducing r||"“ we avoid the region more affected by shot 
noise. As we verified, changing also and r™", or adopting 
different scale selection criteria, do not affect significantly the re¬ 
sults. The aim here is to search for optimal regions in this plane 
to possibly get unbiased constraints. As non-linear dynamics im- 
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Fig. 8. The percentage systematic errors on f(z)(Ti{z) for clus¬ 
ters, as a function of redshift and sample selections. All symbols 
are as in Fig.[2 


pact on different scales for different redshifts and biases of the 
tracers, we expect that the best values of r™" and that is 
the ones that minimise systematics, will be different for different 
sample selections. 

The results of this analysis are shown in Fig. |4] with open 
and solid red dots, obtained by htting the data with the model 
given by Eqs. (Unil and (O, respectively, that is by considering 
either the linear or the non-linear RSD model. We explore the 
ranges 5 < Mpc] < 20 and 15 < Mpc] < 35, 

highlighted by the grey areas in the bottom panels of the lower 
window of Fig.|4l Specihcally, we consider a 10x10 grid in the 
r™", plane. The best-ht values of r™" and shown in 
the bottom panels of the lower window, are the ones that min¬ 
imise the systematic error on /erg. The ranges considered are 
large enough to get systematic errors on /erg lower than ~ 5%. 
Indeed, for proper values of FJ*” and it is possible to signifi¬ 
cantly reduce the systematic error on /erg at all redshifts, without 
having to improve the treatment of non-linearities. The error bars 
have been estimated using the scaling formula given by Eq.[T9l 

Both the linear and non-linear dispersion models provide 
similar results when and are kept free. In other words, it 
seems more convenient simply to not consider the scales where 
non-linear effects have a non-negligible impact, rather than to 
try modelling them with the empirical description provided by 
Eq. ([ISll. 

These results show that i) it is possible to reduce signifi¬ 
cantly the systematics in RSD constraints, even with the disper¬ 
sion model, by cutting the plane conveniently, and ii) 

the optimal FJ*”, range depends on sample selection, thus it 
cannot be just fixed in mulfi-fracer analyses. As it can be seen 
in Fig. m we do not hnd any clear trend of Fjj” and as a 
function of redshift, as it would be expected if the systematics 
were caused by non-linearities. This is possibly caused by the 
not sufficiently large volumes considered. The analysis of larger 


Fig. 9. The iso-correlation contours of the redshift- 
space two-point correlation function, corresponding to 
the values /(^±,^||)= [0.1,0.2,0.4,1,3], for AGN with 

\og{MBH[ft ^ Mq]) > 6.3, at z = 0.2 (black contours). The 
other contours show the best-fit models, as in Fig. [3] 

data sets, that we defer to a future work, will hopefully clarify 
this point. 

For completeness, we then apply our method to several 
subsamples with different selections. Specihcally, we consider 
galaxy catalogues selected in stellar mass, Mstar, g-band abso¬ 
lute magnitude. Mg, and star formation rate, SFR. The selection 
thresholds considered and the number of galaxies in each sam¬ 
ple are reported in Table Ol The percentage systematic errors 
on /cTg are reported in Fig.|5] as a function of redshift and sam¬ 
ple selections, as indicated by the labels. The result is quite re¬ 
markable: it is indeed possible to get almost unbiased constraints 
on /(Tg with the dispersion model, independent of sample selec¬ 
tions, provided that r™" and r™' are chosen conveniently. For 
the largest galaxy sample shown in Fig. [2l the preferred value 
of r™" is ~ 15 h^'Mpeat low and high redshifts, and ~ 5 

Mpe around z = 1, while is ~ 15 Mpeat low red¬ 
shifts, and increases up to ~ 35 Mpe at z = 2. Again, we do 
not hnd any clear trend of FJ*“ and r™”* as a function of redshift 
or sample bias. 

Overall, all of these results show that our modelling of RSD 
is mostly sensitive to the lower htting cut-off, depending on how 
the tracers relate to the underlying mass. Indeed, the effect of 
the different selections considered here is just to select subsam¬ 
ples of haloes of different masses (hence bias) hosting the ob¬ 
served galaxies. As noted above, this has to be considered in 
multi-tracer analyses, where a single sample selection might in¬ 
troduce systematics. 

5.2.2. Galaxy groups and clusters 

In this section, we perform the same analysis presented in ^5.2.11 
on the mock samples of galaxy groups and clusters. To simplify 
the discussion, in the following we will use the term clusters to 


9 





















































F. Marulli et al.: Redshift-space distortions with different tracers 



0 0.5 1 1.5 2 

redshift 



0.5 1 1.5 2 

redshift 



Fig. 10. The best-fit values of f(z)crs(z) for AGN with 
log(MBH[f!^* Mq]) > 6.3. All the symbols are as in Fig.|4l 


refer to all of these objects, though the less massive ones should 
be considered as galaxy groups, or just haloes, from an observa¬ 
tional perspective (see Table lC.2]) . 

Fig. 0 shows the iso-correlation contours of the redshift- 
space two-point correlation function of clusters with 
log(M 5 oo[/t^* Mq]) > 13, at z = 0.2. The other contours 
are the best-fit models, as in Fig. [3] As it can be seen, the 
fingers-of-God anisotropies are almost absent, differently from 
the galaxy correlation function shown in Fig.[3j due to the lower 
values of non-linear motions of clusters at small scales. This 
makes the convolution of Eq. (fTbl) negligible, as it can be seen 
comparing the dashed blue and solid red lines. Interestingly, 
when the real-space correlation function ^(r) is measured from 
the mocks (green contours), the dispersion model fails to match 
the small-scale clustering shape. This is primarily caused by 
the paucity of the cluster sample. Due to the low number of 
cluster pairs at small separations, the small-scale clustering 
cannot be estimated accurately, thus introducing systematics in 
the model. On the contrary, when ^(r) is estimated from CAMB, 
the small-scale clustering shape can be accurately described. 
For visual clarity, we have shown here only the iso-correlation 
contours down to f(ij_, iy) = 0.2, as for lower values they appear 
too scattered. 


Fig. 11. The percentage systematic errors on f(z)crs(z) for AGN, 
as a function of redshift and sample selections. All symbols are 
as in Fig. [5] 


Fig. Q shows the best-fit values of /erg as a function of red¬ 
shift. In contrast with what we found with galaxies, the con¬ 
straints on the linear growth rate obtained in the scale ranges 
3 < Mpe] < 35 and 3 < Mpe] < 35 (blue squares) 
result severely overestimated at all redshifts, independent of us¬ 
ing the linear or the non-linear dispersion model. To investigate 
how these systematics depend on the scale range considered in 
the fit, we repeat the same analysis described in ^5.2.11 keeping 
free the parameters and The result is shown by the red 
dots in Fig.|7] As one can see, we are able to substantially reduce 
the systematics on /erg, at all redshifts. This is obtained by cut¬ 
ting the small perpendicular separations, limiting the analysis at 
rj^ > 20 hr^ Mpe, as shown in the lower window of Fig.]?] 

These result s appear in overall ag reement with what previ¬ 
ously found by [Bianchi et al.l (120121) where systematic errors 
become positive for DM halo masses larger than 10'^/;“' Mq, 
that is for group masses and above. This represents a nice con¬ 
firmation, not related to the specific simulation or halo/cluster 
selection, but only to the dynamics of haloes of large masses. 

Then, we consider several cluster subsamples selected in 
mass, M 500 , temperature, Tsoo, a nd lum inosity, L 500 , whose main 
properties are reported in Table ICjR Fig. [8] shows the percent¬ 
age systematic errors on /erg as a function of redshift and sample 
selections. For the largest samples analysed we do not get any 
significant bias on /erg, similarly to what we found with galax¬ 
ies. However, when selecting too few objects, the results appear 
quite scattered, though without systematic trends. Indeed, the 
advantages resulting from the low fingers-of-God anisotropies 
compensate with the paucity of cluster pairs at small scales. 
Summarising, it seems possible to get unbiased constraints on 


® M 500 . T 500 and L 500 are, respectively, the mass, the temperature and 
the luminosity enclosed within a sphere of radius rjooc, in which the 
mean matter density is equal to 500 times the critical density pcrit(z) - 
3 N^(z)I87 !'G, where ff(z) is the Huhble parameter. 
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the linear growth rate from the RSD of galaxy clusters, but the 
cluster sample has to be sufficiently numerous. 

5.2.3. AGN 

Finally, we analyse the AGN mock samples. The results are 
shown in Figs.[9llT0landlTn that are the analogous of the plots 
presented in the previous two sections for galaxies and clus¬ 
ters. As expected, we find similar results to the ones found with 
galaxies. Differently from the cluster case, the fingers-of-God 
are clearly evident in the AGN correlation function, and not ac¬ 
curately described by the dispersion model, as it can be seen in 
Fig. 0 Nevertheless, the convolution of Eq. (fTbl) significantly 
helps in reducing the systematic error on /erg, even more than 
what we found when analysing the galaxy sample, as evident 
by comparing Figs. [3 and 0 Also in this case, it is possible 
to get almost unbiased constraints on fcr% at all redshifts and 
for the different selections considered, provided that the anal¬ 
ysis is restricted at FJ*” > 15 hT^ Mpe and < 15 Mpe. 
Specifically, we considered here AGN catalogues selected in BH 
mass, Mbh, bolometric luminosity, Lboi, and Eddington factor, 
/Edd, as reported in Table Ic3] Only the AGN samples with the 
lower number of objects analysed here (cyan points in Fig. [TTb 
are not large enough to provide robust constraints on /erg. 


6. Discussion 

The results presente d in this paper extend the previous work by 
[Bianchi et akl (l2012h . who analysed mock samples of DM haloes 
at z = 1. One of the main differences with respect to that work 
is that here we used galaxy, cluster and AGN simulated cata¬ 
logues, that are more directly comparable to observational sam¬ 
ples. Moreover, we investigated the systematic errors of /erg, 
instead of yS, and extended the analysis at different redshifts and 
sample selections. The use of these simulated objects as test par¬ 
ticles provides a better treatment of non-linearities in small-scale 
dynamics, that translates to a more realistic description of the 
fingers-of-God in the redshift-space galaxy and AGN clustering, 
as it can be seen in Figs.[3]and|9] On the other hand, the selected 
cluster s amples are more sim ilar to the DM halo samples anal¬ 
ysed by [Bianchi et al.l (l2012l) . as expected, with fingers-of-God 
anisotropies almost absent, as shown in Fig.jb] Nevertheless, the 
Magneticum samples provide cluster properties that allowed us to 
apply more realistic selections. 

One of the main results found by [Bianchi et akl (12012h is 
that, when using the dispersion model, the systematic errors on 
the linear growth rate are positive for DM haloes with masses 
lower than ~ 10*^/; * Mp, and negative o therwise (see also 
lOkumura & Jindl201 ll : iMarulli et al. l2012bli . Our findings re¬ 
inforce these conclusions. Indeed, regardless of the specific se¬ 
lections considered, the overall trends of our systematic errors 
depend ultimately on the mass of the DM haloes hosting the se¬ 
lected tracers. 

Moreover, extending the [Bianchi et al.1 (12012h analysis, we 
investigated how the systematic errors depend on the comov¬ 
ing scales considered. Our main finding is that it is possible to 
substantially reduce the systematic errors at all redshifts and se¬ 
lections by restricting the analysis in a proper subregion of the 
r± - r\\ plane. The latter depends on the properties of the selected 
objects, specifically on the mass of the host halo, due to the dif¬ 
ficulties in modelling the non-Gaussian nature of the velocity 
probability density function of the tracers. 


7. Conclusions 

Clustering anisotropies in redshift-space are one of the most 
effective probes to test Einstein’s General Relativity on large 
scales. Though only galaxy samples have been considered so 
far for these analyses, other astronomical tracers, such as clus¬ 
ters and AGN, will be used in the next future, aimed at max¬ 
imising the dynamic and redshift ranges where to constrain the 
linear growth rate of cosmic structures. In particular, clusters 
of galaxies can be efficiently used to probe the low redshift 
Universe. Their large bias and low velocity dispersion at small 
scal es make them optimal tracers for large scale structur e analy¬ 
ses (ISereno et al.ll2015l:IVeror)alumbo et^l20l4l2016ll . On the 
other hand, AGN samples can be exploited to push the GR test 
on larger redshifts. 

In this work, we made use of both galaxy, cluster and AGN 
mock samples extracted from the Magneticum, a cosmological 
hydrodinamic simulation of the standard ACDM Universe, to 
investigate the accuracy of the widely used dispersion model as 
a function of scale and sample selection. Instead of analysing 
the clustering multipoles, we preferred to extract /erg constraints 
from the monopole only, or from the full 2D anisotropic correla¬ 
tion /(ij., i'll). This allowed us to explore the rj^ - ry plane, aimed 
at finding the optimal range of scales to minimise systematic un¬ 
certainties. 

The main results of this analysis are the following. 

- It is possible to get almost unbiased constraints on fcr% by 
modelling the large-scale monopole of the two-point corre¬ 
lation function with the Kaiser limit of the dispersion model 
(Eq. I 1 Ibl) . provided that the linear bias factor bcr% is known. 
The latter has to be estimated at sufficiently large scales, to 
avoid non-linearities. With real data, the linear bias can be 
estimated either with the deprojection technique, or from the 
three-point correlation function, or by combining clustering 
and lensing measurements (see the discussion in ® . 

- When higher multipoles are considered, /erg and bcr% can 
be jointly constrained, and the s tatistical error on ferg 
is substantially red uced (see e.g. Ide la Torre et al.l 1201 3t 
[Petracca et al.ll201^ . Flowever, systematic errors can arise 
if non-linearities are not treated accurately, such as with the 
dispersion model (Eas. fTOlfTTl l. On the other hand, the anal¬ 
ysis shown in this paper demonstrates that it is still possible 
to reduce the systematic errors on fcr% if the range of scales 
used to model RSD is chosen appropriately. The latter de¬ 
pends on the properties of the tracers. We do not find how¬ 
ever any clear prescription to relate these scale ranges to the 
sample properties. 

- If the linear growth rate is estimated from the RSD of either 
galaxies or AGN, the value of /erg will be underestimated at 
z < 1, and overestimated at larger redshifts, when the fit is 
performed in the fixed range of scales 3 < rj^, r\\[h^^ Mpe] < 
35. With rf" ~ 15 r'Mpe and rj'” ~ 20 r'Mpc, we 
get almost unbiased constraints. However, the proper val¬ 
ues of PJ*™ and depend on sample selection, and should 
be fixed differently for each sample analysed, calibrating the 
analysis using specific mock catalogues. 

- If the RSD of galaxy groups and clusters are modelled in 

the range of scales 3 < r\\{hr^ Mpe] < 35, the value of 

/cTg results severely overestimated at all redshifts. The re- 
su lts at z = 1 appear in overall agreement with what found 
bv iBianchi et al.1 (120121) . The systematic error can be reduced 
if the fit is performed at larger scales, as expected due to the 
low number of cluster pairs at small separations. Specifically, 


11 





































F. Marulli et al.: Redshift-space distortions with different tracers 


with r™" ~ 20 Mpc, we get unbiased constraints, but 

only for the densest group samples. The paucity of the most 

massive cluster samples does not allow us to obtain robust 

results. 

To summarise, as the dispersion model cannot accurately 
describe the non-linear motions at small scales, the constraints 
on /cTg can be biased. At the present time, it is difficult to in¬ 
clude model uncertainties in this kind of analyses, and even 
more to improve the RSD modelling at small scales (see e.g. 
Ide la Torre & Guz^ 120121 : iBianchi et alJ 120151 and references 
therein). Instead of trying to improve the model, in this analysis 
we investigated how the systematic errors on the linear growth 
rate caused by model uncertainties depend on the comoving 
scales considered. As we have shown, in order to reduce these 
systematics it is necessary to restrict the analysis in a proper sub- 
region of the rj_ - ry plane. More specifically, we found that it is 
enough to choose a proper value for r™" and r™“. Indeed, as we 
verified, changing also and r|j"“ does not affect signihcantly 
the results. 

Overall, we hnd that, in order to reduce the systematics in the 
/cTg measurements, it is more convenient simply to not consider 
the small scales where non-linear effects distort significantly the 
clustering shape, rather than to try modelling them with the em¬ 
pirical description provided by the dispersion model given by 
Eq. (fTbll . The drawback of this approach is to inevitably increase 
the statistic error on /erg. Moreover, the - ry region to be se¬ 
lected depends on the properties of the tracers, and should be 
determined using suitable mock catalogues. On the other hand, 
without a sufficiently accurate modelisation of non-linear dy¬ 
namics at small scales, this approach is the only one possible 
to get unbiased constraints. 

This can hav e a non-negligible impact on multiple tracer 
analyses (see e.g. [McDonald & Seli^l2009l : iBlake et ^120131: 

iMarfn et'^l2016l and references therein) if the RSD are mod¬ 

elled at small non-linear scales with the dispersion model. To 
avoid different systematic errors on the linear growth rate esti¬ 
mated from the small-scale RSD of different tracers, the rj_ - ry 
plane analysed should be different for each sample. On the other 
hand, in multiple tracer techniques that require to analyse differ¬ 
ent samples at the same scales, it is necessary to limit the analy¬ 
sis at large enough scales where systematic errors are lower than 
statistical errors for all the different tracers. And the same when 
joining together different redshifts, where RSD constraints can 
be differently biased. 
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Fig. 12. The number of mock galaxies {red dots), clusters {blue 
squares) and AGN {green diamonds) as a function of redshift, 
for the three selections: Gsl, Cll, Afl (see Tables ICTl 1C. 21 and 
01 ) . 
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Appendix A: The main properties of the seiected 
samples 

Tables [cm IC. 21 and IC.3l rer)ort the main properties of the mock 
samples analysed in this work. We consider only integral cata¬ 
logues. The quantities min and max reported in the Tables are 
the minimum and maximum thresholds used for the selections, 
while median is the median of the sample distributions. As it 
can be seen, we adopt suitable sample selections to have a min¬ 
imum of ~ 50000 objects in the smallest galaxy samples, and 
~ 10000 objects in the smallest cluster and AGN samples. As an 
illustrative example. Fig. [12] shows the number of mock galax¬ 
ies, clusters and AGN as a function of redshift, for three different 
selections, as indicated by the labels. 

Appendix B: The effect of random catalogue 
number densities 

We test the impact of our assumptions in constructing the ran¬ 
dom catalogues by repeating the full statistical analysis on a sub¬ 
set of the mock samples, using random catalogues with different 
number densities. Fig. lB.ll shows an example case obtained from 
the sample of mock galaxies at z = 0.2. The lines show the ratio 
A^/cr^ as a function of scales, where = ^Nr=z - with Nr 
being the ratio between the number of random objects and galax¬ 
ies, and cr^ is the estimated statistical uncertainty. The different 
lines refer to the values = 1,5,10. The effect of different 
random number densities on the measured correlation functions 
is marginal, considering the estimated uncertainties. As verified, 
the net effect on the final RSD constraints is also negligible. 
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20 40 60 80 100 

r 

Fig.B.l. The ratio Af/cr^ as a function of scales obtained from 
the sample of mock galaxies at z = 0.2, where = ftvR =3 “ 
with Nr being the ratio between the number of random objects 
and galaxies, and cr^ is the estimated statistical uncertainty. Blue, 
red and green lines refer to the values of = 1,5,10, respec¬ 
tively. 


0 0.5 1 1.5 2 

redshift 

Fig. C.l. The difference between fcr% measurements obtained 
with diagonal Poisson errors and by considering the full covari¬ 
ance matrix assessed with either the jackknife (solid lines) or 
the bootstrap (dashed lines) techniques, divided by the Icr un¬ 
certainty on /cTg, as a function of redshift. Red, blue and green 
symbols refer to galaxy, cluster and AGN samples, respectively. 


Appendix C: The effect of covariance matrix 
approximations 

All the results presented in this paper have been obtained by as¬ 
suming Poisson errors in clustering measurements and by ne¬ 
glecting the off-diagonal elements of the covariance matrix. Both 
of these approximations can introduce systematics in the inferred 
cosmological constraints. However, as we verified, they do not 
impact significantly our present conclusions, introducing just 
marginal differences in the final obtained constraints, consider¬ 
ing the estimated uncertainties. Fig. 1C.II shows the difference 
in the fcr% measurements obtained with diagonal Poisson errors 
and by considering the full covariance matrix assessed with ei¬ 
ther jackknife or bootstrap techniques, divided by the Icr uncer¬ 
tainty on /cTg. We use 125 iackknife/bootstr ap mock realisations 
to co mpute the covariance matrices (see IVeropalumbo et ahl 
1201 6l for more details on the used methods and codes). As it 
can be seen, the effect of assuming negligible off-diagonal el¬ 
ements in the covariance matrix is not statistically significant. 
Moreover, the off-diagonal elements in the covariance matrix 
can introduce systematics and spurious scatter in the inferred 
cosmological constraints, if not properly estimated with a suffi¬ 
ciently large number of mocks. The above considerations moti¬ 
vated our choice of neglecting off-diagonal terms in the analyses 
pre sented in this paper, in line with previous similar works (see 
e.g. lBianchi et al.ll201.5l i lPetracca et al.ll201fih . 
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Table C.l. Main properties of the samples of galaxies analysed in the paper. 


redshift log(AfsTAR[/i ‘ Mp]) Mg log(5Ffi[Mo/yegr]) 



sample 

min 

max 

median 

objects 

sample 

min 

max 

median 

objects 

sample 

min 

max 

median 

objects 

0.2 

Gml 

10.00 

13.24 

10.78 

2565412 

Gmgl 

-20.00 

-27.39 

-21.90 

2555977 

Gsl 

1.00 

3.40 

1.38 

1460374 


Gm2 

10.45 


10.98 

1881270 

Gmg2 

-20.94 


-21.93 

2465212 

Gs2 

1.24 


1.47 

1057208 


Gm3 

10.90 


11.13 

1080529 

Gmg3 

-21.87 


-22.31 

1333968 

Gs3 

1.48 


1.64 

509094 


Gm4 

11.35 


11.55 

255117 

Gmg4 

-22.81 


-23.19 

267639 

Gs4 

1.72 


1.86 

169018 


Gm5 

11.79 


11.96 

50000 

Gmg5 

-23.74 


-24.12 

49995 

Gs5 

1.96 


2.08 

49996 

0.52 

Gml 

10.00 

13.10 

10.73 

2167365 

Gmgl 

-20.00 


-21.65 

2154603 

Gsl 

1.00 

3.33 

1.45 

1432387 


Gm2 

10.42 


10.93 

1556875 

Gmg2 

-20.88 


-21.70 

1985640 

Gs2 

1.27 


1.54 



Gm3 

10.84 


11.11 

905634 

Gmg3 

-21.76 


-22.21 

913480 

Gs3 

1.54 


1.71 

550535 


Gm4 

11.27 


11.46 

245678 

Gmg4 

-22.64 



231625 

Gs4 

1.80 


1.95 

181393 


Gm5 

11.69 


11.84 

50000 

Gmg5 

-23.52 


-23.88 

49993 

Gs5 



2.20 

49998 

0.72 

Gml 

10.00 

13.03 

10.69 

1894509 

Gmgl 

-20.00 

-26.78 

-21.53 

1872485 

Gsl 

1.00 

3.35 

1.50 

1360951 


Gm2 

10.40 


10.90 

1347360 

Gmg2 

-20.86 


-21.62 

1657323 

Gs2 

1.28 


1.58 

1069332 


Gm3 

10.81 


11.10 

788388 

Gmg3 

-21.72 


-22.19 

708654 

Gs3 

1.56 


1.75 

569550 


Gm4 

11.21 


11.39 

242502 

Gmg4 

-22.58 



209604 

Gs4 

1.84 


2.00 

190113 


Gm5 

11.61 


11.77 

50000 

Gmg5 

-23.44 


-23.79 

49984 

Gs5 

2.13 


2.25 

49999 

1 

Gml 

10.00 

12.92 

10.62 

1507943 

Gmgl 

-20.00 

-26.78 

-21.43 

1465209 

Gsl 

1.00 

3.41 

1.58 

1206017 


Gm2 

10.37 


10.84 


Gmg2 



-21.55 

1229380 

Gs2 

1.30 


1.65 

994050 


Gm3 

10.75 


11.07 

620010 

Gmg3 

-21.69 


-22.15 

490496 

Gs3 

1.60 


1.81 

572970 


Gm4 

11.12 


11.28 

247210 

Gmg4 

-22.53 


-23.09 

171265 

Gs4 

1.89 


2.05 

199961 


Gm5 

11.49 


11.65 

49999 

Gmg5 

-23.38 


-23.71 

49689 

Gs5 

2.19 


2.32 

50000 

1.5 

Gml 

10.00 

12.79 

10.55 

1028296 

Gmgl 

-20.00 


-21.28 

959724 

Gsl 

1.00 

3.53 

1.67 

921472 


Gm2 

10.33 


10.76 

711625 

Gmg2 

-20.79 


-21.44 

760311 

Gs2 

1.31 


1.72 

811432 


Gm3 

10.67 


10.99 

421555 

Gmg3 

-21.58 


-21.96 

298389 

Gs3 

1.62 


1.87 

516807 


Gm4 

11.00 


11.19 

204649 

Gmg4 

-22.38 


-23.28 

86610 

Gs4 

1.93 


2.09 

206667 


Gm5 

11.33 


11.48 

50000 

Gmg5 

-23.17 


-23.56 

49994 

Gs5 

2.24 


2.37 

50000 

2 

Gml 

10.00 

12.49 

10.48 

618406 

Gmgl 

-20.00 

-25.87 

-21.17 

552844 

Gsl 

1.00 

3.54 

1.74 

589434 


Gm2 

10.29 


10.67 

425026 

Gmg2 

-20.51 


-21.27 

477671 

Gs2 

1.31 


1.78 

543876 


Gm3 

10.58 


10.88 

254679 

Gmg3 



-21.48 


Gs3 

1.61 


1.90 

386208 


Gm4 

10.87 


11.10 


Gmg4 

-21.53 


-21.85 

148563 

Gs4 

1.92 


2.11 

179827 


Gm5 

11.16 



49999 

Gmg5 

-22.05 


-22.49 

47889 

Gs5 

2.23 


2.36 

49996 
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Table C.2. Main properties of the samples of galaxy groups and clusters analysed in the paper. 


redshift XogjMsmlh ' Mp]) TsmikeV] log(L 5 oo[erg/j])~ 



sample 

min 

max 

median 

objects 

sample 

min 

max 

median 

objects 

sample 

min 

max 

median 

objects 

0.2 

Cml 

13.00 

14.97 

13.26 

103201 

Ctl 

0.10 

8.27 

0.29 

274787 

Cll 

42.00 

46.57 

42.51 

249655 


Cm2 

13.18 


13.41 

63705 

Ct2 

0.19 


0.35 

215729 

C12 

42.41 


42.75 

153119 


Cm3 

13.37 


13.58 

36858 

Ct3 

0.34 


0.54 

110534 

C13 

42.83 


43.16 

65040 


Cm4 

13.55 


13.74 

19963 

Ct4 

0.64 


0.89 

39922 

C14 

43.24 


43.54 

26922 


Cm5 

13.74 


13.90 

9999 

Ct5 

1.18 


1.52 

9902 

C15 

43.66 


43.92 

9977 

0.52 

Cml 

13.00 

14.97 

13.23 

84523 

Ctl 

0.10 

7.94 

0.30 

247341 

Cll 

42.00 

46.09 

42.64 

239836 


Cm2 

13.16 


13.37 

53798 

Ct2 

0.18 


0.36 

198050 

C12 

42.43 


42.78 

180659 


Cm3 

13.31 


13.50 

32449 

Ct3 

0.33 


0.52 

109150 

C13 

42.87 


43.19 

75378 


Cm4 

13.47 


13.64 

18518 

Ct4 

0.61 


0.84 

40267 

C14 

43.30 


43.59 

29210 


Cm5 

13.62 


13.79 

10000 

Ct5 

1.11 


1.42 

9959 

C15 

43.73 


43.98 

9991 

0.72 

Cml 

13.00 

14.69 

13.22 

69404 

Ctl 

0.10 

6.57 

0.31 

218241 

Cll 

42.00 

46.24 

42.77 

215861 


Cm2 

13.13 


13.33 

45804 

Ct2 

0.18 


0.36 

176278 

C12 

42.44 


42.83 

187288 


Cm3 

13.27 


13.45 

29018 

03 

0.32 


0.50 

103154 

C13 

42.89 


43.20 

84083 


Cm4 

13.40 


13.57 

17408 

Ct4 

0.58 


0.80 

39329 

C14 

43.33 


43.61 

31049 


Cm5 

13.54 


13.69 

9999 

05 

1.04 


1.31 

9973 

C15 

43.77 


44.02 

9992 

1 

Cml 

12.50 

14.58 

12.81 

167773 

Ctl 

0.10 

5.85 

0.31 

166270 

Cll 

42.00 

46.15 

42.88 

166445 


Cm2 

12.72 


12.97 

104160 

Ct2 

0.17 


0.36 

137909 

C12 

42.44 


42.91 

156759 


Cm3 

12.95 


13.15 

55569 

Ct3 

0.30 


0.47 

86893 

C13 

42.88 


43.20 

84056 


Cm4 

13.17 


13.34 

25461 

Ct4 

0.52 


0.72 

36186 

C14 

43.32 


43.59 

32064 


Cm5 

13.39 


13.54 

10000 

Ct5 

0.91 


1.15 

9991 

C15 

43.76 


43.99 

9992 

1.5 

Cml 

12.50 

14.43 

12.77 

100290 

Ctl 

0.10 

4.83 

0.31 

99747 

Cll 

42.00 

46.20 

43.06 

100220 


Cm2 

12.67 


12.90 

66549 

Ct2 

0.16 


0.35 

86504 

C12 

42.44 


43.06 

98730 


Cm3 

12.85 


13.04 

39631 

Ct3 

0.27 


0.45 

57585 

C13 

42.87 


43.23 

70071 


Cm4 

13.02 


13.18 

21126 

Ct4 

0.44 


0.61 

30005 

C14 

43.31 


43.60 

29699 


Cm5 

13.20 


13.34 

10000 

Ct5 

0.72 


0.92 

9975 

C15 

43.75 


43.97 

9971 

2 

Cml 

12.50 

14.23 

12.74 

46857 

Ctl 

0.10 

4.17 

0.28 

46602 

Cll 

42.00 

45.79 

43.21 

46857 


Cm2 

12.61 


12.82 

34568 

Ct2 

0.15 


0.31 

42502 

C12 

42.40 


43.21 

46818 


Cm3 

12.73 


12.91 

24316 

Ct3 

0.22 


0.39 

30646 

C13 

42.80 


43.23 

44871 


Cm4 

12.84 


13.00 

15934 

Ct4 

0.34 


0.52 

18884 

C14 

43.20 


43.52 

23962 


Cm5 

12.95 


13.10 

10000 

Ct5 

0.50 


0.68 

9934 

C15 

43.60 


43.87 

9953 
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Table C.3. Main properties of the samples of AGN analysed in the paper. 


redshift log(MBH[/i ‘ Mp]) log(Lboi[erg/j]) log(/Edd) 



sample 

min 

max 

median 

objects 

sample 

min 

max 

median 

objects 

sample 

min 

max 

median 

objects 

0.2 

Ami 

6.30 

11.05 

8.50 

946365 

All 

40.00 

48.21 

43.26 

903502 

Afl 

-5.00 

0.19 

-3.28 

900571 


Am2 

7.16 


8.76 

559616 

A12 

41.49 


43.35 

820496 

Af2 

-3.78 


-2.82 

669586 


Am3 

8.02 


8.78 

530118 

A13 

42.99 


43.67 

561712 

Af3 

-2.57 


-1.65 

282446 


Am4 

8.88 


9.13 

200602 

A14 

44.48 


45.26 

104962 

Af4 

-1.35 


-0.92 

100757 


Am5 

9.74 


9.88 

9999 

A15 

45.97 


46.29 

9999 

Af5 

-0.13 


0.11 

9999 

0.52 

Ami 

6.30 

11.03 

8.49 

687068 

All 

40.00 

48.27 

43.55 

662433 

Afl 

-5.00 

0.20 

-2.98 

664774 


Am2 

7.15 


8.84 

387939 

A12 

41.55 


43.65 

602637 

Af2 

-3.74 


-2.75 

559834 


Am3 

8.00 


8.86 

368515 

A13 

43.10 


43.93 

435109 

Af3 

-2.47 


-1.57 

227878 


Am4 

8.85 


9.12 

189324 

A14 

44.64 


45.39 

86203 

Af4 

-1.21 


-0.69 

74218 


Am5 

9.70 


9.83 

10000 

A15 

46.19 


46.51 

9999 

Af5 

0.05 


0.18 

9999 

0.72 

Ami 

6.30 

10.98 

8.18 

556411 

All 

40.00 

48.37 

43.65 

539859 

Afl 

-5.00 

0.20 

-2.80 

542532 


Am2 

7.14 


8.89 

298180 

A12 

41.54 


43.76 

493217 

Af2 

-3.75 


-2.64 

483058 


Am3 

7.99 


8.92 

282007 

A13 

43.09 


44.07 

351857 

Af3 

-2.50 


-1.63 

216373 


Am4 

8.83 


9.13 

170321 

A14 

44.63 


45.37 

85212 

Af4 

-1.25 


-0.70 

69112 


Am5 

9.67 


9.80 

9999 

A15 

46.18 


46.44 

9999 

Af5 

0.00 


0.16 

10000 

1 

Ami 

6.30 

10.92 

7.99 

362513 

All 

40.00 

48.40 

43.91 

353674 

Afl 

-5.00 

0.19 

-2.63 

355647 


Am2 

7.13 


9.02 

194441 

A12 

41.60 


44.03 

323269 

Af2 

-3.72 


-2.50 

324605 


Am3 

7.97 


9.05 

181542 

A13 

43.20 


44.36 

235508 

Af3 

-2.43 


-1.46 

154064 


Am4 

8.80 


9.17 

138296 

A14 

44.80 


45.56 

69691 

Af4 

-1.15 


-0.46 

56167 


Am5 

9.63 


9.75 

9999 

A15 

46.39 


46.65 

10000 

Af5 

0.14 


0.18 

9959 

1.5 

Ami 

6.30 

10.75 

7.01 

199546 

All 

40.00 

48.52 

44.17 

196211 

Afl 

-5.00 

0.20 

-2.05 

197164 


Am2 

7.12 


9.14 

97867 

A12 

41.68 


44.30 

181622 

Af2 

-3.71 


-1.90 

186168 


Am3 

7.94 


9.19 

89539 

A13 

43.36 


44.70 

133756 

Af3 

-2.41 


-1.05 

114389 


Am4 

8.75 


9.27 

73532 

A14 

45.04 


46.14 

49508 

Af4 

-1.12 


-0.29 

60067 


Am5 

9.57 


9.69 

9999 

A15 

46.72 


46.95 

10000 

Af5 

0.17 


0.18 

8465 

2 

Ami 

6.30 

10.64 

6.68 

88522 

All 

40.00 

48.63 

44.09 

87627 

Afl 

-5.00 

0.20 

-1.55 

87913 


Am2 

7.08 


9.13 

38422 

A12 

41.66 


44.27 

82555 

Af2 

-3.76 


-1.50 

85568 


Am3 

7.86 


9.22 

33150 

A13 

43.31 


45.13 

56194 

Af3 

-2.53 


-1.13 

67129 


Am4 

8.64 


9.32 

27031 

A14 

44.97 


46.31 

30473 

Af4 

-1.30 


-0.40 

37606 


Am5 

9.41 


9.57 

10001 

A15 

46.63 


46.92 

10000 

Af5 

-0.06 


0.18 

9988 


16 



F. Marulli et al.: Redshift-space distortions with diiferent tracers 


References 

Amendola, L. et al. 2013, Living Reviews in Relativity, 16, 6 
Bender, F. et al. 2012, Mon. Not. R. Astron. Soc., 423, 3430 
Bianchi, D., Chiesa, M., & Guzzo, L. 2015, Mon. Not. R. Astron. Soc., 446, 75 
Bianchi, D., Guzzo, L., Branchini, E., et al. 2012, Mon. Not. R. Astron. Soc., 
427, 2420 

Blake, C. et al. 2011, Mon. Not. R. Astron. Soc., 415, 2876 
Blake, C. et al. 2012, Mon. Not. R. Astron. Soc., 425, 405 
Blake, C. et al. 2013, Mon. Not. R. Astron. Soc., 436, 3089 
Bocquet, S., Saro, A., Dolag, K., & Mohr, J. J. 2016, Mon. Not. R. Astron. Soc., 
456, 2361 

Bondi, H. 1952, Mon. Not. R. Astron. Soc., 112, 195 
Bondi, H. & Hoyle, F. 1944, Mon. Not. R. Astron. Soc., 104, 273 
Bonoli, S., Marulli, R, Springel, V., et al. 2009, Mon. Not. R. Astron. Soc., 396, 
423 

Cabre, A. & Gaztanaga, E. 2009a, Mon. Not. R. Astron. Soc., 393, 1183 
Cabre, A. & Gaztanaga, E. 2009b, Mon. Not. R. Astron. Soc., 396, 1119 
Chabrier, G. 2003, Publ. Astron. Soc. Pac., 115, 763 

Chuang, C.-H., Prada, R, Cuesta, A. J., et al. 2013, Mon. Not. R. Astron. Soc., 
433,3559 

Chuang, C.-H. & Wang, Y. 2013, Mon. Not. R. Astron. Soc., 435, 255 
Contreras, C. et al. 2013, Mon. Not. R. Astron. Soc., 430, 924 
Davis, M. & Peebles, P. J. E. 1983, Astrophys. J., 267, 465 
de la Torre, S. & Guzzo, L. 2012, Mon. Not. R. Astron. Soc., 427, 327 
de la Torre, S. et al. 2013, Astron. Astrophys., 557, A54 
Di Matteo, T, Springel, V., & Hemquist, L. 2005, Nature, 433, 604 
Dolag, A. et al. in preparation 

Dolag, K., Vazza, E, Brunetti, G., & Tormen, G. 2005, Mon. Not. R. Astron. 
Soc., 364, 753 

Fabjan, D., Borgani, S., Tomatore, L., et al. 2010, Mon. Not. R. Astron. Soc., 
401, 1670 

Ferland, G. J., Korista, K. T., Vemer, D. A., et al. 1998, Publ. Astron. Soc. Pac., 
no, 761 

Fisher, K. B., Scharf, C. A., & Lahav, O. 1994, Mon. Not. R. Astron. Soc., 266, 
219 

Guzzo, L. et al. 2008, Nature, 451, 541 

Haardt, F. & Madau, P. 2001, in Clusters of Galaxies and the High Redshift 
Universe Observed in X-rays, ed. D. M. Neumann & J. T. V. Tran, 64 
Hamilton, A. J. S. 1992, Astrophys. J. Lett., 385, L5 

Hamilton, A. J. S. 1998, in Astrophysics and Space Science Library, Vol. 231, 
The Evolving Universe, ed. D. Hamilton, 185 
Hawkins, E. et al. 2003, Mon. Not. R. Astron. Soc., 346, 78 
Hearin, A. P. 2015, Mon. Not. R. Astron. Soc., 451, L45 
Hirschmann, M., Dolag, K., Saro, A., et al. 2014, Mon. Not. R. Astron. Soc., 
442, 2304 

Hoyle, F. & Lyttleton, R. A. 1939, Proceedings of the Cambridge Philosophical 
Society, 35, 405 

Kaiser, N. 1987, Mon. Not. R. Astron. Soc., 227, 1 
Komatsu, E. et al. 2011, Astrophys. J. Suppl., 192, 18 
Landy, S. D. & Szalay, A. S. 1993, Astrophys. J., 412, 64 
Laureijs, R. et al. 2011, ArXiv e-prints, astro-ph/1110.3193 
Lewis, A. & Bridle, S. 2002, Phys. Rev. D, 66, 103511 
Lilje, P. B. & Efstathiou, G. 1989, Mon. Not. R. Astron. Soc., 236, 851 
Macaulay, E., Wehus, I. K., & Eriksen, H. K. 2013, Physical Review Letters, 
111, 161301 

Marin, F. A., Beutler, E, Blake, C., et al. 2016, Mon. Not. R. Astron. Soc., 455, 
4046 

Marulli, E, Baldi, M., & Moscardini, L. 2012a, Mon. Not. R. Astron. Soc., 420, 
2377 

Marulli, E, Bianchi, D., Branchini, E., et al. 2012b, Mon. Not. R. Astron. Soc., 
426, 2566 

Marulli, E, Bonoli, S., Branchini, E., et al. 2009, Mon. Not. R. Astron. Soc., 396, 
1404 

Marulli, F., Bonoli, S., Branchini, E., Moscardini, L., & Springel, V. 2008, Mon. 
Not. R. Astron. Soc., 385, 1846 

Marulli, E, Carbone, C., Viel, M., Moscardini, L., & Cimatti, A. 2011, Mon. 
Not. R. Astron. Soc., 418, 346 

Marulli, R, Veropalumbo, A., & Moresco, M. 2016, Astronomy and Computing, 
14, 35 

McDonald, P. & Seljak, U. 2009, J. Cosm. Astro-Particle Phys., 10, 7 
McGill, C. 1990, Mon. Not. R. Astron. Soc., 242, 428 
Merloni, A. et al. 2012, ArXiv e-prints, astro-ph/1209.3114 
Mo, H. J., Jing, Y. P, & Boemer, G. 1992, Astrophys. J., 392, 452 
Mohammad, F. G., de la Torre, S., Bianchi, D., Guzzo, L., & Peacock, J. A. 2016, 
Mon. Not. R. Astron. Soc., 458, 1948 

Moresco, M., Marulli, R, Baldi, M., Moscardini, L., & Cimatti, A. 2014, Mon. 
Not. R. Astron. Soc., 443, 2874 


Okumura, T. & Jing, Y. P. 2011, Astrophys. J., 726, 5 

Padovani, P. & Matteucci, F. 1993, Astrophys. J., 416, 26 

Peacock, J. A. & Dodds, S. J. 1996, Mon. Not. R. Astron. Soc., 280, L19 

Peacock, J. A. et al. 2001, Nature, 410, 169 

Peebles, P. J. E. 1980, The large-scale structure of the universe 

Percival, W. J. et al. 2004, Mon. Not. R. Astron. Soc., 353, 1201 

Perlmutter, S. et al. 1999, Astrophys. J., 517, 565 

Petracca, E, Marulli, E, Moscardini, L., et al. 2016, Mon. Not. R. Astron. Soc., 
462, 4208 

Reid, B. A., Seo, H.-J., Leauthaud, A., Tinker, J. L., & White, M. 2014, Mon. 
Not. R. Astron. Soc., 444, 476 

Reid, B. A. et al. 2012, Mon. Not. R. Astron. Soc., 426, 2719 
Riess, A. G. et al. 1998, Astron. J., 116, 1009 
Ross, N. P. et al. 2007, Mon. Not. R. Astron. Soc., 381, 573 
Samushia, L., Percival, W. J., & Raccanelli, A. 2012, Mon. Not. R. Astron. Soc., 
420, 2102 

Saro, A. et al. 2014, Mon. Not. R. Astron. Soc., 440, 2610 
Scoccimarro, R. 2004, Phys. Rev. D, 70, 083007 
Seljak, U. & McDonald, P. 2011, J. Cosm. Astro-Particle Phys., 11, 39 
Sereno, M., Veropalumbo, A., Marulli, E, et al. 2015, Mon. Not. R. Astron. Soc., 
449, 4147 

Shakura, N. I. & Sunyaev, R. A. 1973, Astron. Astrophys., 24, 337 
Smith, R. E., Peacock, J. A., Jenkins, A., et al. 2003, Mon. Not. R. Astron. Soc., 
341, 1311 

Spergel, D. et al. 2013, ArXiv e-prints, astro-ph/1305.5422 
Springel, V. 2005, Mon. Not. R. Astron. Soc., 364, 1105 
Springel, V. & Hemquist, L. 2002, Mon. Not. R. Astron. Soc., 333, 649 
Springel, V. & Hemquist, L. 2003a, Mon. Not. R. Astron. Soc., 339, 289 
Springel, V. & Hemquist, L. 2003b, Mon. Not. R. Astron. Soc., 339, 312 
Springel, V. et al. 2005, Nature, 435, 629 

Taruya, A., Nishimichi, T., & Saito, S. 2010, Phys. Rev. D, 82, 063522 
Tegmark, M. et al. 2004, Astrophys. J., 606, 702 

Thielemann, F.-K., Argast, D., Brachwitz, F., et al. 2003, Nuclear Physics A, 
718, 139 

Tojeiro, R. et al. 2012, Mon. Not. R. Astron. Soc., 424, 2339 
Tomatore, L., Borgani, S., Dolag, K., & Matteucci, F. 2007, Mon. Not. R. Astron. 
Soc., 382, 1050 

van den Hoek, L. B. & Groenewegen, M. A. T. 1997, Astron. Astrophys. Suppl., 
123, 305 

Veropalumbo, A., Marulli, E, Moscardini, L., Moresco, M., & Cimatti, A. 2014, 
Mon. Not. R. Astron. Soc., 442, 3275 

Veropalumbo, A., Marulli, E, Moscardini, L., Moresco, M., & Cimatti, A. 2016, 
Mon. Not. R. Astron. Soc., 458, 1909 

Wang, L., Reid, B., & White, M. 2014, Mon. Not. R. Astron. Soc., 437, 588 
Wiersma, R. P. C., Schaye, J., & Smith, B. D. 2009, Mon. Not. R. Astron. Soc., 
393, 99 

Woosley, S. E. & Weaver, T. A. 1995, Astrophys. J. Suppl., 101, 181 
Zhang, H., Yu, H., Noh, H., & Zhu, Z.-H. 2008, Physics Letters B, 665, 319 
Zurek, W. H., Quinn, P. J., Salmon, J. K., & Warren, M. S. 1994, Astrophys. J., 
431, 559 

Zwicky, F. 1937, Astrophys. J., 86, 217 


17 



